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We have performed muon spin rotation/relaxation (/.tSR) measurements on single crystals of the 
chiral helimagnet Cr 1 / 3 NbS 2 at zero to low magnetic field. The transition from the paramagnetic 
to helical magnetically ordered phase at zero held is marked by the onset of a coherent oscillation 
of the zero-held muon spin polarization below a critical temperature T c . An enhancement of the 
muon spin precession frequency is observed below T ~ 50 K, where anomalous behavior has been 
observed in bulk transport measurements. The enhanced precession frequency indicates a low- 
temperature modification of the helical magnetic structure. A Landau free energy analysis suggests 
that the low-temperature change in the magnetic structure is caused by a structural change. We also 
suggest a longer periodicity of helicity below T~ 50 K, which can be verihed by neutron scattering 
experiments. 

PACS numbers: 75.30.-m, 75.10.-b, 76.75.+i 


Cr!/ 3 NbS 2 is an itinerant chiral helimagnet 
(CHM) with an hexagonal non-centrosymmetric 
crystal structured In crystals belonging to a non- 
centrosymmetric space group, CHM order is generally 
understood to be a consequence of competing symmetric 
Heisenberg exchange and Dzyaloshinsky-Moriya (DM) 
interactions. The CHM order in Cr 1 / 3 NbS 2 continuously 
evolves into an incommensurate chiral magnetic solition 
lattice (CSL) upon application of a magnetic field 
perpendicular to the helical c-axis (H_Lc)dd The CSL 
consists of periodic domains of ferromagnetically ordered 
spins in the ab plane separated by 360° spin domain 
walls. The period of the CSL increases monotonically 
with the magnitude of the applied field, resulting 
in a continuous phase transition to a commensurate 
ferromagnetic (FM) state. The ability to easily tune 
the size of the magnetic domains of the CSL, and hence 
the magnetic potential experienced by the spins of the 
intinerant electrons, makes Cr 1 / 3 NbS 2 an appealing 
candidate for spintronics applications. 

At elevated temperatures there is a transition to a 
paramagnetic (PM) state. Reported values of the PM- 
to-CHM phase transition temperature T c vary between 
118 K (Ref. 0 ) and 133.5 K (Ref. 0 ). While the long- 
range magnetically ordered phases of Cr!/ 3 NbS 2 have 
been established, the details of the transition from the 
PM phase to the CHM or CSL phases are unresolved. 
Anomalies in the magnetic field dependence of the resis¬ 
tivity, and the temperature dependences of the Seebeck 
and Hall coefficients near T c have been reported.® - ' In 
particular, at T c a signficant negative magnetoresistance 
(with Hie) has been observed, with no saturation of the 
resistivity at fields where both the resistivity and mag¬ 
netization saturate at lower temperatures^ The Seebeck 
and Hall coefficients also show a pronounced minimum 


near T c , indicative of changes to the electronic structure 
and/or the effects of the magnetism on the mobility of 
the charge carriers. Likewise, anomalous behavior has 
been observed in magnetoresistance, thermopower and 
Hall resistivity data below T~ 50 K£ At these low tem¬ 
peratures the magnetoresistance at high fields becomes 
anisotropic, there is a rapid decrease of the Seebeck coef¬ 
ficient with decreasing temperature, and the field depen¬ 
dence of the Hall resistivity drastically changes. To gain 
futher insight into the sources of the transport anoma¬ 
lies, we have used the /rSR technique to investigate the 
nature of the accompanying magnetism. 

Our /iSR measurements were performed at TRIUMF 
on a ~ 50 mm 2 mosaic consisting of seven ~ 0.25 mm 
thick plate-like Cr 1 / 3 NbS 2 single crystals, which were 
grown as described in Ref. @. The positive muon (/r + ) 
beam was directed at the large face of the mosaic, with 
the linear momentum of the beam parallel to the c-axis 
of the crystals (z-direction). The initial muon spin po¬ 
larization P(f = 0) for all measurements was rotated to 
be perpendicular to the c-axis (x-direction). Data was 
collected with zero applied magnetic field (ZF), and for 
a transverse field (TF) [i.e. perpendicular to P(f = 0)] 
applied perpendicular to the c-axis (y-direction). Note 
that for the desired Hlc orientation the field must be 
applied perpendicular to the beam, since the crystals are 
too thin to orient with the beam shining on the a-c or b-c 
plane faces. The ensuing beam deflection by the magnetic 
force exerted on the incoming /.t + limited the applied field 
strength to H< 600 Oe. 

Figure [IJ a) shows representative ZF-/rSR asymmetry 
spectra. The data for 135 <T < 170 K were best fit to 
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FIG. 2: (Color online) Simulations of the sample (s) contribu¬ 
tion to the ZF-/rSR signal for Cr!/ 3 NbS 2 in the CHM ground 
state. The dashed curves show the contribution from muons 
stopping near Nb at (0.77, 0.32, 0.96) and near S at (0.37, 
0.78, 0.6), where the coordinate values are multiples of the 
respective lattice constants a = b= 5.7 A, and c=12.1 A. The 
Cr magnetic moment is assumed to be 4.4 /r_g. The occupa¬ 
tion probability of the muon site near S is taken to be twice 
that of the site near Nb. 


where G'kt(A,, t.) is a static Gaussian Kubo-Toyabe re¬ 
laxation function, and Ab is a small time-independent 
background signal. The fits yield A\/A 2 ~ 2, Ai = 
0.091(4) /us -1 , A 2 = 0.414(5) /us -1 , and a relaxation rate 
A that exhibits a divergent increase as the temperature is 
lowered. Based on these results and nuclear dipole field 
calculations, we attribute the two sample components to 
a muon stopping site near a S atom and a site near Nb. 
The amplitude ratio A 1 /A 2 implies that the probability 
of the muon stopping at the S site is approximately twice 
that of the Nb site. 

For T < 130 K the relaxation rate is sufficiently fast 
that some of the measured ZF-/uSR signal is damped out 
in the initial dead time of the /iSR spectromter, resulting 
in a reduction of the initial amplitude Aq = A\ + A 2 + Ab 
[see Fig.[ljc)]. The asymmetry spectra at these temper¬ 
atures were fit to 


A(t) = Aie Alt cos( 27 t ft) + A 2 e As< + Ab . (2) 


FIG. 1: (Color online) (a) Representative ZF-pSR asymmetry 
spectra of Cr 1/ / 3 NbS 2 recorded at different temperatures. The 
inset shows the asymmetry spectra at early times, with the 
data packed into smaller time bins. The solid curves through 
the data points are fits to Eqs. © or ©. (b) Representa¬ 
tive TF-/rSR asymmetry spectra of Cr!/ 3 NbS 2 recorded at 
different temperatures with H = 600 Oe. The solid curves 
through the data points are fits to Eq. ©. (c) Temperature 
dependence of the total initial asymmetry Ao. 


The oscillating term is indicative of magnetic order, and 
the spin precession frequency / is related to the magni¬ 
tude of the mean field B^ sensed by the muon ensemble, 
where / = (y M /2tt)B^ and 7 M /27r = 135.54 MHz/T is the 
muon gyromagnetic ratio. Despite there being two muon 
stopping sites, we could not resolve two distinct frequen¬ 
cies. This is understandable from simulations of the con¬ 
tribution of each of these sites to the ZF-/rSR signal (see 
Fig- O, where the signal from muons stopping near Nb 
is smaller and rapidly damped. The temperature depen¬ 
dence of / indicates a magnetic transition at T c = 129.6 K 
(see Fig. [3]) , which is close to previous estimates of the 
Curie temperature of Cr 1 / 3 NbS 2 >i 
Figure ©b) shows representative TF-/iSR spectra for 
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H _L c and H = 600 Oe. The external field adds vectorially 
to the local field generated by the spin structure. At all 
temperatures the TF-p,SR spectrum is well described by 
the sum of sample and background contribution of the 
following form 

A(t) = A s e~ Kt cos( 27 t ft) + A^e -17 * cos(27r/bf), (3) 

where Ao=A s + Ab is the total signal amplitude. There 
is a rapid signal depolarization below T c , which mani¬ 
fests itself in a reduction of Aq that is more severe than 
with H = 0 [see Fig. [TJc)] . Bulk magnetization mea¬ 
surements suggest that the CSL state is not fully formed 
until H « 850 Oe.— A distortion of the CHM order into an 
imhomogeneous helicoid is consistent with the enhanced 
depolarization observed here at weaker fields. As shown 
in Fig.[3](a), the temperature dependence of the spin pre¬ 
cession frequency / is similar to that observed with H = 0, 
but shifted in value by the applied field. There is some 
enhancement of / at low temperatures, which is indica¬ 
tive of a structural change and/or modification of the 
spin structure. At H = 0 Oe and H = 100 Oe this occurs 
below T « 50 K, where as mentioned earlier, anomalous 
behavior has been observed in transport measurements. 

Figure (3][b) shows the temperature dependence of the 
frequency difference f(T) — f(T = 170 K). With H > 0, 
the spin precession frequency increases as the zero-field 
magnetic transition temperature T c is approached from 
above. As shown in the inset of Fig. [3j[b) , f(T) — f(T = 
170 K) oc H at temperatures T c < T < 170 K, where 
the proportionality constant is roughly proportional to 
1/(T — T c ). Hence the observed frequency shift above T c 
is consistent with the low-field response of a paramagnet, 
whereby the induced magnetiztion is linearly dependent 
on H. A full and temperature-independent amplitude Aq 
is observed above T c [see Fig. [TJb)], indicating that the 
paramagnetism occurs throughout the sample volume. 

To provide a possible explanation for the enhanced 
mean magnetic field sensed by the muon ensemble (~/) 
at low temperatures, we adopt a Landau free energy 
analysis. In the absence of the H _L c magnetic field, 
the helimagnetic phase occurs via FM and DM interac¬ 
tions due to a lack of inversion symmetry. As shown in 
Fig. U the Cr atoms form layered triangular lattices in 
Cr!/ 3 NbS 2 , with an hexagonal space group of P6 3 22. For 
each Cr atom in a layer, there are three nearest-neighbor 
Cr atoms in an adjacent layer. C 2 rotational symme¬ 
try exists about an axis located halfway between two Cr 
atoms in the different layers, with inversion symmetry 
about this point broken by the S atoms. The DM vector 
linking the two Cr atoms (Di) then lies in a plane per¬ 
pendicular to the C 2 axis, as shown in Fig. [4j Since the 
three Cr atoms within a layer are related by C 3 rotational 
symmetry (about an axis parallel to the 2 -axis), there 
are three DM vectors related by symmetry (Di, D 2 , and 
D 3 ), which when summed together give a resultant DM 
vector D pointing along the 2 -axis. Consequently, the 
spin density of the helical order can be described by&2 

S(r) = S exp(tQ • r) + S* exp(—iQ ■ r), (4) 




FIG. 3: (Color online) Temperature dependence of (a) the 
muon spin precession frequency / observed for different ap¬ 
plied magnetic field strengths, and of (b) the frequency shift 
f(T) — f(T =170 K) plotted with a logarithmic vertical scale. 
The black solid curves through the data points in both pan¬ 
els are guides to the eye. The red curve in (b), which ex¬ 
tends off the logarithmic scale, comes from a fit of the high- 
temperature data for H = 0 Oe to f(T) = /(0)(1 — T/T c ) n , 
yielding T c = 129.58(4) K and n = 0.366(6). The inset shows 
the magnetic field dependence of f(T)—f(T =170 K) at tem¬ 
peratures above T c . 

where Q = Qqz , S = S x x+iS y y, and by sysmmetry | S x \ = 
|5 , y | = S'. Since the free energy should be invariant under 
the exchange of Q —» Q and S —> — S*, it is given to 
fourth order by^ 

F = a 0 S ■ S* + DQ 0 S 2 + cnQlS 2 + PS 4 , (5) 

where D represents the strength of the DM interaction^ 
and the parameter /3 is assumed to be positive, such that 
the free energy is bound. Minimizing F with respect to 
Q 0 , we find Qo = —D/2ai, resulting in the simplified free 
energy equation 

F = a Q S 2 + pS 4 , (6) 

where «q =oto — D 2 /4a\. Subsequent minimization of 
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FIG. 4: (Color online) Illustration of the DM interaction in 
Crj/ 3 NbS 2 . Each Cr atom (blue spheres) is linked to three Cr 
atoms in one of the adjacent planes (indicated by the blue ar¬ 
rows). The S and Nb atoms are the yellow and green spheres, 
respectivley. The cylindrical red arrow represents a C 2 ro¬ 
tational symmetry axis located halfway between the two Cr 
atoms. The DM vector Di linking these two Cr atoms (flat 
red arrow) lies in the red-shaded plane perpendicular to the 
C 2 axis. Note the precise orientation of Di in the red-shaded 
plane is unknown. The related DM vectors D 2 and D 3 link¬ 
ing the lower Cr atom with the other two Cr atoms in the 
upper layer are not shown. The gray arrow represents a C 3 
rotational symmetry axis relating the three Cr atoms in the 
upper layer. The short blue arrow above this is the resultant 
DM vector D = Di+D 2 +D 3 . 


S , which implies an increase in the negative value of cxq 
and/or a decrease in the positive value of /3. Since /3 
originates from the spin-spin interaction, smaller /3 im¬ 
plies a stronger FM interaction between the Cr spins - 
although this effect would be quite small. On the other 
hand, a larger negative value of cxq corresponds to a 
weaker DM interaction D. This will occur if there is 
a change in the ratio between the a (and b ) and c-axis 
lattice constants, such that there is a rotation of the vec¬ 
tors Di, D 2 , and D 3 towards the xy plane. Thus we 
propose a change in the length of one or more of the 
lattice constants below T~ 50 K. 

In conclusion, our /iSR measurements reveal a possi¬ 
ble change in the magnetic structure of Cr 1 ( / 3 NbS 2 below 
T ~ 50 K, which is the same temperature region where 
anomalous behavior has been observed in bulk transport 
measurements. The low-temperature modification of the 
magnetic structure can be explained by a change in lat¬ 
tice constants, which can be confirmed by precision mea¬ 
surements of the crystal structure. Moreover, to confirm 
that the change in lattice constants is such that it causes 
a reduction of the DM interaction, we suggest a neutron 
scattering measurement to look for a change of the or¬ 
dering wavevector Q 0 — since a consequence of a smaller 
DM interaction is a longer spin-helix period. Lastly, we 
do not find anything unusual about the magnetism near 
the PM transition at low field (such as phase separation) 
that would contribute to the anomalies oberved in trans¬ 
port measurements near T c . 
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